function y = func_s(p, CONSTS)
    eps_a = CONSTS.eps_a;
    y = (p.^2 - eps_a).^0.5; %sqrt_from_negative_num(p.^2 - eps_a);
end
